Severe adult hemophagocytic lymphohistiocytosis (HLHa) correlates with HLH-related gene variants.

BACKGROUND: The contribution of genetic factors to the severity of adult hemophagocytic lymphohistiocytosis (HLHa) remains unclear. OBJECTIVE: We sought to assess a potential link between HLHa outcomes and HLH-related gene variants. METHODS: Clinical characteristics of 130 HLHa patients (age ≥ 18 years and HScore ≥ 169) and genotype of 8 HLH-related genes (LYST, PRF1, UNC13-D, STX11, STXBP2, RAB27A, XIAP, and SAP) were collected. A total of 34 variants found in only 6 genes were selected on the basis of their frequency and criteria predicted to impair protein function. Severity was defined by refractory disease to HLH treatment, death, or transfer to an intensive care unit. RESULTS: HLHa-associated diseases (ADs) were neoplasia (n = 49 [37.7%]), autoimmune/inflammatory disease (n = 33 [25.4%]), or idiopathic when no AD was identified (n = 48 [36.9%]). Infectious events occurred in 76 (58.5%) patients and were equally distributed in all ADs. Severe and refractory HLHa were observed in 80 (61.5%) and 64 (49.2%) patients, respectively. HScore, age, sex ratio, AD, and infectious events showed no significant association with HLHa severity. Variants were identified in 71 alleles and were present in 56 (43.1%) patients. They were distributed as follows: 44 (34.4%), 9 (6.9%), and 3 (2.3%) patients carrying 1, 2, and 3 variant alleles, respectively. In a logistic regression model, only the number of variants was significantly associated with HLHa severity (1 vs 0: 3.86 [1.73-9.14], P = .0008; 2-3 vs 0: 29.4 [3.62-3810], P = .0002) and refractoriness (1 vs 0: 2.47 [1.17-5.34], P = .018; 2-3 vs 0: 13.2 [2.91-126.8], P = .0003). CONCLUSIONS: HLH-related gene variants may be key components to the severity and refractoriness of HLHa.

Allele-dependent interaction of LRRK2 and NOD2 in leprosy.

Leprosy, caused by Mycobacterium leprae, rarely affects children younger than 5 years. Here, we studied a multiplex leprosy family that included monozygotic twins aged 22 months suffering from paucibacillary leprosy. Whole genome sequencing identified three amino acid mutations previously associated with Crohn's disease and Parkinson's disease as candidate variants for early onset leprosy: LRRK2 N551K, R1398H and NOD2 R702W. In genome-edited macrophages, we demonstrated that cells expressing the LRRK2 mutations displayed reduced apoptosis activity following mycobacterial challenge independently of NOD2. However, employing co-immunoprecipitation and confocal microscopy we showed that LRRK2 and NOD2 proteins interacted in RAW cells and monocyte-derived macrophages, and that this interaction was substantially reduced for the NOD2 R702W mutation. Moreover, we observed a joint effect of LRRK2 and NOD2 variants on Bacillus Calmette-Guérin (BCG)-induced respiratory burst, NF-κB activation and cytokine/chemokine secretion with a strong impact for the genotypes found in the twins consistent with a role of the identified mutations in the development of early onset leprosy.

Deep resequencing identifies candidate functional genes in leprosy GWAS loci.

Leprosy is the second most prevalent mycobacterial disease globally. Despite the existence of an effective therapy, leprosy incidence has consistently remained above 200,000 cases per year since 2010. Numerous host genetic factors have been identified for leprosy that contribute to the persistently high case numbers. In the past decade, genetic epidemiology approaches, including genome-wide association studies (GWAS), identified more than 30 loci contributing to leprosy susceptibility. However, GWAS loci commonly encompass multiple genes, which poses a challenge to define causal candidates for each locus. To address this problem, we hypothesized that genes contributing to leprosy susceptibility differ in their frequencies of rare protein-altering variants between cases and controls. Using deep resequencing we assessed protein-coding variants for 34 genes located in GWAS or linkage loci in 555 Vietnamese leprosy cases and 500 healthy controls. We observed 234 nonsynonymous mutations in the targeted genes. A significant depletion of protein-altering variants was detected for the IL18R1 and BCL10 genes in leprosy cases. The IL18R1 gene is clustered with IL18RAP and IL1RL1 in the leprosy GWAS locus on chromosome 2q12.1. Moreover, in a recent GWAS we identified an HLA-independent signal of association with leprosy on chromosome 6p21. Here, we report amino acid changes in the CDSN and PSORS1C2 genes depleted in leprosy cases, indicating them as candidate genes in the chromosome 6p21 locus. Our results show that deep resequencing can identify leprosy candidate susceptibility genes that had been missed by classic linkage and association approaches.

Controlling for human population stratification in rare variant association studies.

Population stratification is a confounder of genetic association studies. In analyses of rare variants, corrections based on principal components (PCs) and linear mixed models (LMMs) yield conflicting conclusions. Studies evaluating these approaches generally focused on limited types of structure and large sample sizes. We investigated the properties of several correction methods through a large simulation study using real exome data, and several within- and between-continent stratification scenarios. We considered different sample sizes, with situations including as few as 50 cases, to account for the analysis of rare disorders. Large samples showed that accounting for stratification was more difficult with a continental than with a worldwide structure. When considering a sample of 50 cases, an inflation of type-I-errors was observed with PCs for small numbers of controls (≤ 100), and with LMMs for large numbers of controls (≥ 1000). We also tested a novel local permutation method (LocPerm), which maintained a correct type-I-error in all situations. Powers were equivalent for all approaches pointing out that the key issue is to properly control type-I-errors. Finally, we found that power of analyses including small numbers of cases can be increased, by adding a large panel of external controls, provided an appropriate stratification correction was used.

Taking population stratification into account by local permutations in rare-variant association studies on small samples.

Many methods for rare variant association studies require permutations to assess the significance of tests. Standard permutations assume that all individuals are exchangeable and do not take population stratification (PS), a known confounding factor in genetic studies, into account. We propose a novel strategy, LocPerm, in which individual phenotypes are permuted only with their closest ancestry-based neighbors. We performed a simulation study, focusing on small samples, to evaluate and compare LocPerm with standard permutations and classical adjustment on first principal components. Under the null hypothesis, LocPerm was the only method providing an acceptable type I error, regardless of sample size and level of stratification. The power of LocPerm was similar to that of standard permutation in the absence of PS, and remained stable in different PS scenarios. We conclude that LocPerm is a method of choice for taking PS and/or small sample size into account in rare variant association studies.

Genome-wide association study of resistance to Mycobacterium tuberculosis infection identifies a locus at 10q26.2 in three distinct populations.

The natural history of tuberculosis (TB) is characterized by a large inter-individual outcome variability after exposure to Mycobacterium tuberculosis. Specifically, some highly exposed individuals remain resistant to M. tuberculosis infection, as inferred by tuberculin skin test (TST) or interferon-gamma release assays (IGRAs). We performed a genome-wide association study of resistance to M. tuberculosis infection in an endemic region of Southern Vietnam. We enrolled household contacts (HHC) of pulmonary TB cases and compared subjects who were negative for both TST and IGRA (n = 185) with infected individuals (n = 353) who were either positive for both TST and IGRA or had a diagnosis of TB. We found a genome-wide significant locus on chromosome 10q26.2 with a cluster of variants associated with strong protection against M. tuberculosis infection (OR = 0.42, 95%CI 0.35-0.49, P = 3.71×10-8, for the genotyped variant rs17155120). The locus was replicated in a French multi-ethnic HHC cohort and a familial admixed cohort from a hyper-endemic area of South Africa, with an overall OR for rs17155120 estimated at 0.50 (95%CI 0.45-0.55, P = 1.26×10-9). The variants are located in intronic regions and upstream of C10orf90, a tumor suppressor gene which encodes an ubiquitin ligase activating the transcription factor p53. In silico analysis showed that the protective alleles were associated with a decreased expression in monocytes of the nearby gene ADAM12 which could lead to an enhanced response of Th17 lymphocytes. Our results reveal a novel locus controlling resistance to M. tuberculosis infection across different populations.

The complex pattern of genetic associations of leprosy with HLA class I and class II alleles can be reduced to four amino acid positions.

Leprosy is a chronic disease caused by Mycobacterium leprae. Worldwide, more than 200,000 new patients are affected by leprosy annually, making it the second most common mycobacterial disease after tuberculosis. The MHC/HLA region has been consistently identified as carrying major leprosy susceptibility variants in different populations at times with inconsistent results. To establish the unambiguous molecular identity of classical HLA class I and class II leprosy susceptibility factors, we applied next-generation sequencing to genotype with high-resolution 11 HLA class I and class II genes in 1,155 individuals from a Vietnamese leprosy case-control sample. HLA alleles belonging to an extended haplotype from HLA-A to HLA-DPB1 were associated with risk to leprosy. This susceptibility signal could be reduced to the HLA-DRB1*10:01~ HLA-DQA1*01:05 alleles which were in complete linkage disequilibrium (LD). In addition, haplotypes containing HLA-DRB3~ HLA-DRB1*12:02 and HLA-C*07:06~ HLA-B*44:03~ HLA-DRB1*07:01 alleles were found as two independent protective factors for leprosy. Moreover, we replicated the previously associated HLA-DRB1*15:01 as leprosy risk factor and HLA-DRB1*04:05~HLA-DQA1*03:03 as protective alleles. When we narrowed the analysis to the single amino acid level, we found that the associations of the HLA alleles were largely captured by four independent amino acids at HLA-DRß1 positions 57 (D) and 13 (F), HLA-B position 63 (E) and HLA-A position 19 (K). Hence, analyses at the amino acid level circumvented the ambiguity caused by strong LD of leprosy susceptibility HLA alleles and identified four distinct leprosy susceptibility factors.

Family-based genome-wide association study of leprosy in Vietnam.

Leprosy is a chronic infectious disease of the skin and peripheral nerves with a strong genetic predisposition. Recent genome-wide approaches have identified numerous common variants associated with leprosy, almost all in the Chinese population. We conducted the first family-based genome-wide association study of leprosy in 622 affected offspring from Vietnam, followed by replication in an independent sample of 1181 leprosy cases and 668 controls of the same ethnic origin. The most significant results were observed within the HLA region, in which six SNPs displayed genome-wide significant associations, all of which were replicated in the independent case/control sample. We investigated the signal in the HLA region in more detail, by conducting a multivariate analysis on the case/control sample of 319 GWAS-suggestive HLA hits for which evidence for replication was obtained. We identified three independently associated SNPs, two located in the HLA class I region (rs1265048: OR = 0.69 [0.58-0.80], combined p-value = 5.53x10-11; and rs114598080: OR = 1.47 [1.46-1.48], combined p-value = 8.77x10-13), and one located in the HLA class II region (rs3187964 (OR = 1.67 [1.55-1.80], combined p-value = 8.35x10-16). We also validated two previously identified risk factors for leprosy: the missense variant rs3764147 in the LACC1 gene (OR = 1.52 [1.41-1.63], combined p-value = 5.06x10-14), and the intergenic variant rs6871626 located close to the IL12B gene (OR = 0.73 [0.61-0.84], combined p-value = 6.44x10-8). These results shed new light on the genetic control of leprosy, by dissecting the influence of HLA SNPs, and validating the independent role of two additional variants in a large Vietnamese sample.

Genome-wide association study of Buruli ulcer in rural Benin highlights role of two LncRNAs and the autophagy pathway.

Buruli ulcer, caused by Mycobacterium ulcerans and characterized by devastating necrotizing skin lesions, is the third mycobacterial disease worldwide. The role of host genetics in susceptibility to Buruli ulcer has long been suggested. We conduct the first genome-wide association study of Buruli ulcer on a sample of 1524 well characterized patients and controls from rural Benin. Two-stage analyses identify two variants located within LncRNA genes: rs9814705 in ENSG00000240095.1 (P = 2.85 × 10-7; odds ratio = 1.80 [1.43-2.27]), and rs76647377 in LINC01622 (P = 9.85 × 10-8; hazard ratio = 0.41 [0.28-0.60]). Furthermore, we replicate the protective effect of allele G of a missense variant located in ATG16L1, previously shown to decrease bacterial autophagy (rs2241880, P = 0.003; odds ratio = 0.31 [0.14-0.68]). Our results suggest LncRNAs and the autophagy pathway as critical factors in the development of Buruli ulcer.

Skin-specific antibodies neutralizing mycolactone toxin during the spontaneous healing of Mycobacterium ulcerans infection.

Buruli ulcer, a neglected tropical infectious disease, is caused by Mycobacterium ulcerans. Without treatment, its lesions can progress to chronic skin ulcers, but spontaneous healing is observed in 5% of cases, suggesting the possible establishment of a host strategy counteracting the effects of M. ulcerans. We reveal here a skin-specific local humoral signature of the spontaneous healing process, associated with a rise in antibody-producing cells and specific recognition of mycolactone by the mouse IgG2a immunoglobulin subclass. We demonstrate the production of skin-specific antibodies neutralizing the immunomodulatory activity of the mycolactone toxin, and confirm the role of human host machinery in triggering effective local immune responses by the detection of anti-mycolactone antibodies in patients with Buruli ulcer. Our findings pave the way for substantial advances in both the diagnosis and treatment of Buruli ulcer in accordance with the most recent challenges issued by the World Health Organization.

A 1-Year Prospective French Nationwide Study of Emergency Hospital Admissions in Children and Adults with Primary Immunodeficiency.

PURPOSE: Patients with primary immunodeficiency (PID) are at risk of serious complications. However, data on the incidence and causes of emergency hospital admissions are scarce. The primary objective of the present study was to describe emergency hospital admissions among patients with PID, with a view to identifying "at-risk" patient profiles. METHODS: We performed a prospective observational 12-month multicenter study in France via the CEREDIH network of regional PID reference centers from November 2010 to October 2011. All patients with PIDs requiring emergency hospital admission were included. RESULTS: A total of 200 admissions concerned 137 patients (73 adults and 64 children, 53% of whom had antibody deficiencies). Thirty admissions were reported for 16 hematopoietic stem cell transplantation recipients. When considering the 170 admissions of non-transplant patients, 149 (85%) were related to acute infections (respiratory tract infections and gastrointestinal tract infections in 72 (36%) and 34 (17%) of cases, respectively). Seventy-seven percent of the admissions occurred during winter or spring (December to May). The in-hospital mortality rate was 8.8% (12 patients); death was related to a severe infection in 11 cases (8%) and Epstein-Barr virus-induced lymphoma in 1 case. Patients with a central venous catheter (n = 19, 13.9%) were significantly more hospitalized for an infection (94.7%) than for a non-infectious reason (5.3%) (p = 0.04). CONCLUSION: Our data showed that the annual incidence of emergency hospital admission among patients with PID is 3.4%. The leading cause of emergency hospital admission was an acute infection, and having a central venous catheter was associated with a significantly greater risk of admission for an infectious episode.

Pleiotropic effects for Parkin and LRRK2 in leprosy type-1 reactions and Parkinson's disease.

Type-1 reactions (T1R) are pathological inflammatory episodes and main contributors to nerve damage in leprosy. Here, we evaluate the genewise enrichment of rare protein-altering variants in 7 genes where common variants were previously associated with T1R. We selected 474 Vietnamese leprosy patients of which 237 were T1R-affected and 237 were T1R-free matched controls. Genewise enrichment of nonsynonymous variants was tested with both kernel-based (sequence kernel association test [SKAT]) and burden methods. Of the 7 genes tested 2 showed statistical evidence of association with T1R. For the LRRK2 gene an enrichment of nonsynonymous variants was observed in T1R-free controls (PSKAT-O = 1.6 × 10-4). This genewise association was driven almost entirely by the gain-of-function variant R1628P (P = 0.004; odds ratio = 0.29). The second genewise association was found for the Parkin coding gene PRKN (formerly PARK2) where 7 rare variants were enriched in T1R-affected cases (PSKAT-O = 7.4 × 10-5). Mutations in both PRKN and LRRK2 are known causes of Parkinson's disease (PD). Hence, we evaluated to what extent such rare amino acid changes observed in T1R are shared with PD. We observed that amino acids in Parkin targeted by nonsynonymous T1R-risk mutations were also enriched for mutations implicated in PD (P = 1.5 × 10-4). Hence, neuroinflammation in PD and peripheral nerve damage due to inflammation in T1R share overlapping genetic control of pathogenicity.

French-style genetics v. 2.0: The "e-CohortE" project.

In the digital age, a genetics cohort has become much more than a simple means of determining the cause of a disease. Two-sided markets, of which 23andMe, Ancestry DNA and MyHeritage are the best known, have showed this perfectly over the last few years: a cohort has become a means of producing massive amounts of data for medical, scientific and commercial exploitation, and for genetic use in particular. French law does not currently allow these foreign private companies to develop on French national territory and also forbids the creation of similar entities in France. However, at least in theory, this same law does not preclude the creation of new types of cohorts in France inspired by the success of two-sided markets but retaining features specific to the French healthcare management system. We propose an optimal solution for France, for genomic studies associated with multi-subject questionnaires, still purely theoretical for the moment: the development, with no need for any change in the law, of France's own version of "Genetics v.2.0": "e-CohortE."

Microdeletion on chromosome 8p23.1 in a familial form of severe Buruli ulcer.

Buruli ulcer (BU), the third most frequent mycobacteriosis worldwide, is a neglected tropical disease caused by Mycobacterium ulcerans. We report the clinical description and extensive genetic analysis of a consanguineous family from Benin comprising two cases of unusually severe non-ulcerative BU. The index case was the most severe of over 2,000 BU cases treated at the Centre de Dépistage et de Traitement de la Lèpre et de l'Ulcère de Buruli, Pobe, Benin, since its opening in 2003. The infection spread to all limbs with PCR-confirmed skin, bone and joint infections. Genome-wide linkage analysis of seven family members was performed and whole-exome sequencing of both patients was obtained. A 37 kilobases homozygous deletion confirmed by targeted resequencing and located within a linkage region on chromosome 8 was identified in both patients but was absent from unaffected siblings. We further assessed the presence of this deletion on genotyping data from 803 independent local individuals (402 BU cases and 401 BU-free controls). Two BU cases were predicted to be homozygous carriers while none was identified in the control group. The deleted region is located close to a cluster of beta-defensin coding genes and contains a long non-coding (linc) RNA gene previously shown to display highest expression values in the skin. This first report of a microdeletion co-segregating with severe BU in a large family supports the view of a key role of human genetics in the natural history of the disease.

An eQTL variant of ZXDC is associated with IFN-γ production following Mycobacterium tuberculosis antigen-specific stimulation.

There is a large inter-individual variability in the response to Mycobacterium tuberculosis infection. In previous linkage analyses, we identified a major locus on chromosome region 8q controlling IFN-γ production after stimulation with live BCG (Bacillus Calmette-Guérin), and a second locus on chromosome region 3q affecting IFN-γ production triggered by the 6-kDa early secretory antigen target (ESAT-6), taking into account the IFN-γ production induced by BCG (IFNγ-ESAT6BCG). High-density genotyping and imputation identified ~100,000 variants within each linkage region, which we tested for association with the corresponding IFN-γ phenotype in families from a tuberculosis household contact study in France. Significant associations were replicated in a South African familial sample. The most convincing association observed was that between the IFNγ-ESAT6BCG phenotype and rs9828868 on chromosome 3q (p = 9.8 × 10-6 in the French sample). This variant made a significant contribution to the linkage signal (p < 0.001), and a trend towards the same association was observed in the South African sample. This variant was reported to be an eQTL of the ZXDC gene, biologically linked to monocyte IL-12 production through CCL2/MCP1. The identification of rs9828868 as a genetic driver of IFNγ production in response to mycobacterial antigens provides new insights into human anti-tuberculosis immunity.

Age-Dependent Association of TNFSF15/TNFSF8 Variants and Leprosy Type 1 Reaction.

A current major challenge in leprosy control is the prevention of permanent disabilities. Host pathological inflammatory responses termed type 1 reaction (T1R) are a leading cause of nerve damage for leprosy patients. The environmental or inherited factors that predispose leprosy cases to undergo T1R are not known. However, studies have shown an important contribution of host genetics for susceptibility to T1R. We have previously identified variants encompassing the TNFSF15/TNFSF8 genes as T1R risk factors in a Vietnamese sample and replicated this association in a Brazilian sample. However, we failed to validate in Brazilian patients the strong association of TNFSF15/TNFSF8 markers rs6478108 and rs7863183 with T1R that we had observed in Vietnamese patients. Here, we investigated if the lack of validation of these variants was due to age-dependent effects on association using four independent population samples, two from Brazil and two from Vietnam. In the combined analysis across the four samples, we observed a strong association of the TNFSF15/TNFSF8 variants rs6478108, rs7863183, and rs3181348 with T1R (pcombined = 1.5E-05, pcombined = 1.8E-05, and pcombined = 6.5E-06, respectively). However, the association of rs6478108 with T1R was more pronounced in leprosy cases under 30 years of age compared to the global sample [odds ratio (OR) = 1.95, 95% confidence interval (CI) = 1.54-2.46, pcombined = 2.5E-08 versus OR = 1.46, 95% CI = 1.23-1.73, pcombined = 1.5E-05]. A multivariable analysis indicated that the association of rs6478108 with T1R was independent of either rs7863183 or rs3181348. These three variants are known regulators of the TNFSF8 gene transcription level in multiple tissues. The age dependency of association of rs6478108 and T1R suggests that the genetic control of gene expression varies across the human life span.